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Abstract. Circumstellar discs of Be stars are thought to be formed from material ejected from 
a fast-spinning central star. This material possesses large amounts of angular momentum and 
settles in a quasi-Keplerian orbit around the star. This simple description outlines the basic 
issues that a successful disc theory must address: 1) What is the mechanism responsible for the 
mass ejection? 2) What is the final configuration of the material? 3) How the disc grows? With 
the very high angular resolution that can be achieved with modern interferometers operating in 
Oh the optical and infrared we can now resolve the photosphere and immediate vicinity of nearby 
Be stars. Those observations are able to provide very stringent tests for our ideas about the 
physical processes operating in those objects. This paper discusses the basic hydrodynamics 
of viscous decretion discs around Be stars. The model predictions are quantitatively compared 
to observations, demonstrating that the viscous decretion scenario is currently the most viable 
theory to explain the discs around Be stars. 

£ 

y-n ■ 1. Introduction 

Oh' Recent years witnessed an important progress in our understanding of the circumstellar 
discs of Be stars, largely due to interferometric observations capable of angularly resolve 
those objects at the milliarcsecond (mas) level (see Stee 2010, for a review). 

In the late nineties, discs around Be stars were considered to be equatorially enhanced 
outflowing winds, and several models and mechanisms to drive the outflow were pro- 
posed (see Bjorkman 2000, for a review). The much stronger observational constraints 
t-H ■ available today allow us to rule out several theoretical scenarios that were proposed in 
the past (e.g. the wind compressed disc models of Bjorkman & Cassinelli 1993). As an 
example, spectrointerfometry and spectroastrometry have directly probed the disc kine- 
matics (Meilland et al. 2007b, Stefl et al. 2010, Oudmaijer et al. 2010), revealing that Be 
■ discs rotate very close to Keplerian. As a result, the viscous decretion scenario, proposed 



originally by Lee et al. (1991) and further developed by Porter (1999)[ Okazaki (2001) 



[Bjorkman fc Carciofi (2005) " among others, has emerged as the most viable scenario to 



explain the observed properties of Be discs. 



2. Discs Diagnostics 

Before reviewing the basic aspects of the theory of circumstellar discs, it is useful to 
put in perspective what observations tells us about the structure and kinematics of those 
" ' discs. 



2.1. The formation loci of different observables 

One important issue to consider when analyzing observations is that different observables 
probe different regions of the disc; it is therefore useful to be able to make a correspon- 
dence between a given observable (say, the continuum flux level at a given wavelength) 
and the part of the disc whence it comes. 
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Figure 1. The formation loci of different observables. The calculations assume a rapidly rotating 
BlVe star (T e p ff olc = 25 000 K, fi/n cr it = 0.92) surrounded by viscous decretion discs with 
different sizes. The results correspond to a viewing angle of 30°. Left: continuum emission. 
Plotted is the ratio between the observed flux to the maximum flux, F ma , x , as a function of the 
disc outer radius. F max corresponds to the flux of a model with a disc outer radius of f 000 R. 
Right: polarization and line equivalent width. Calculations were carried out with the hdust code 
(Carciofi & Bjorkman 2006). 



To make such a correspondence, we consider a typical Be star of spectral type BlVe and 
calculate the emergent spectrum arising from viscous decretion discs of different outer 
radii (see Sect.[3]for a detailed description of the structure of viscous discs). In Fig. [T]we 
present results for continuum emission, continuum polarization and the equivalent width 
of emission lines. 

Let us discuss first the continuum and line emission. The reason the spectra of Be 
stars show continuum excess or emission lines is that the dense disc material acts as a 
pseudo-photosphere that is much larger than the stellar photosphere. For the pole-on case 
(inclination i = 0), the effective radius of the pseudo-photosphere is the location where 
the vertical optical depth T\(R c ff) = 1. For continuum Hydrogen absorption, the opacity 
depends on the wavelength as k\ oc A 2 (if one neglects the wavelength dependence of 
the absorption gaunt factors) and thus the size of the pseudo-photosphere increases with 
wavelength (see Carciofi & Bjorkman 2006, appendix A, for a derivation of i? c g(A)). 

In Fig. [TJ left panel, the intersection of the different lines (each corresponding to a 
given wavelength) with the horizontal dashed line marks the position in the disc whence 
about 95% of the continuum excess comes. For instance, from the Figure we see that the 
V band excess is formed very close to the star, within about 2 R, whereas the excess at 1 
mm originates from a much larger volume of the disc. As we shall see below, the fact that 
the continuum excess flux at visible wavelengths forms so close to the star makes such 
observations indispensable to study systems with complex temporal evolution because 
the V band continuum emission has the fastest response to disc changes as a result of 
photospheric activity. 

In the right panel of Fig. [I] we show how the line emission and continuum polarization 
grows with radius. For the optically thick Ha and Br7 lines the disc emission only fills 
in the photospheric absoprtion profile when the disc size is about 5 R. Both lines have 
pseudo-photospheres that extend up to about 20 R. For polarization, the results are at 
odds with the common belief that polarization is formed close to the star; it is seen that 
95% of the maximum polarization is only reached when the disc size is about 10 R. 
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2.2. Disc Thickness 

In a nice example of the diagnostics potential of interferometric and polarimetric obser- 



vations combined, Wood, Bjorkman, & Bjorkman (1997) and Quirrenbach et al. (1997) 



showed that the circumstellar disc of £ Tau is geometrically thin (opening angle of 2.5°). 
This result was confirmed by Carciofi et al. (2009) from a more detailed analysis. Other 



studies based on the fraction of Be-shell stars vs. Be stars typically find larger values for 
the opening angle (e.g. 13° from Hanuschik 1996) but this discrepancy is accommodated 
by the fact that the geometrical thickness of Be discs increases with radius (flared disc) 
and different observables probe different disc regions (Fig. [IJ. In any case, there is no 
doubt that the Be discs are flat, geometrically thin structures. 



2.3. Disc density distribution 

Models derived both empirically (e.g., Waters 1986) and theoretically (e.g., Okazaki 2001 
and Bjorkman & Carciofi 2005) usually predict (or assume) a power-law fall off of the 
disc density \p(r) oc r~ n \. Values for n vary widely in the literature, typically in the range 
2 < n < 4. One question that arises is whether this quoted range for the density slope is 
real or not. 

Below we see that the viscous decretion disc model, in its simplest form of an isothermal 
and isolated disc, predicts a slope of n = 3.5 for the disc density. On one hand, the 
inclusion of other physical effects can change the value of the above slope. For instance, 
non-isothermal viscous diffusion results in much more complex density distribution that 
cannot be well represented by a simple power-law, and tidal effects by a close binary 
may make the density slope shallower (Okazaki et al. 2002). On the other hand, it is 
also true that many of quoted values for n in the literature are heavily influenced by the 
assumptions and methods used in the analysis and should be viewed with caution. 

2.4. Disc dynamics 

Since the earliest detections of Be stars, it became clear that a rotating circumstellar disc- 
like material was the most natural explanation for the observed double-peaked profile of 
emission lines, but the problem of how the disc rotates has been an open issue until 
recently. 

The disc rotation law is profoundly linked with the disc formation mechanism; there- 
fore, determining observationally the disc dynamics is of great interest. Typically, one 
can envisage three limiting cases for the radial dependence of the azimuthal component 
of the velocity, depending on the forces acting on the disc material 

!V X ot ( r / R) ~ 1 radiatively driven outflow , 
V*{r/R) magnetically dominated disc , (2.1) 

Vcrit^R/r) 1 / 2 disc driven by viscosity. 

In the first case (radiatively driven outflow) the dominant force on the material is the 
radially directed radiation pressure that does not exert torques and thus conserves angu- 
lar momentum (V I0 t is thus the rotation velocity of the material when it left the stellar 
surface). The second case corresponds to a rigidly rotating magnestosphere a la a Ori E 
(Townsend et al. 2005), in which the plasma is forced to rotate at the same speed as the 
magnetic field lines (14). The last case corresponds to Keplerian orbital rotation, written 

1/2 

in terms of the critical velocity, V CT it = (GM/R) , which is the Keplerian orbital speed 
at the stellar surface. This case requires a fine-tunning mechanism such that the centrifu- 
gal force, v\jr, exactly balances gravity at all radii. As we shall see below, viscosity does 
provide the fine-tunning mechanism capable of producing a Keplerian disc. 
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Porter & Rivinius (2003) reviewed the then existing observational constraints on the 
disc kinematics and concluded that "all of the kinematic evidence seems to point to 
a disc velocity field dominated by rotation, with little or no radial flow, at least in the 
regions where the kinematic signatures of emission and absorption are significant" . Today, 
spectrointerferometry and spectroastrometry provides clear-cut evidence that, in most 
systems observed and analysed so far (re CMa being the only possible exception, Meilland 
et al. 2007b), the discs rotate in a Keplerian fashion (Meilland et al. 2007b, Stefl et al. 
2010, Oudmaijer et al. 2010). This is an important result, seeing that it indicates that 
viscosity is the driving mechanism of the outflow. 



2.5. Cyclic V/R variations 

About 2/3 of the Be stars present the so-called V/R variations, a phenomenon char- 
acterised by the quasi-cyclic variation in the ratio between the violet and red emission 
peaks of the H I emission lines. These variations are generally explained by global oscilla- 
tions in the circumstellar disc forming a one-armed spiral density pattern that precesses 
around the star with a period of a few years (Kato 1983, Okazaki 1991, Okazaki 1997). 

Recently, Carciofi et al. (2009) provided a quantitative verification of the global disc 
oscillation theory from a detailed modelling of high-angular resolution amber data of 
the Be star £ Tau. From a theoretical perspective, the existence of density waves in 
rotating discs, as suggested by |Kato (1983)| imposes most stringent constraints on the 
rotation velocity, since mode confinement requires that the rotation law be Keplerian 
within about 1% and the radial flow be hightly subsonic (< 0.01 c S) Okazaki 2007). 



2.6. Long-term variations 

One of the most intriguing types of variability observed in the Be stars is the aperiodic 
transition between a normal B phase (discless phase) and a Be phase whereby the disc is 
lost and rebuilt in timescales of months to years (e.g., Clark et al. 2003, Stefl et al. 2010). 
The varying amount of circumstellar gas manifests itself as changes in line profiles (Clark 
et al. 2003), continuum brightness and colors (Harmanec 1983) and polarization (Draper 
et al. 2010). A fine example of the secular process of disc formation and dissipation, 
and its effects on the continuum brightness and colors, is shown in Fig. (5) the outburst 
phase responsible for the disc build-up lasted about 300 hundred days (phases I and II, 
during which the star got brighter and redder) and was followed by a quiescent phase of 
about 500 days (phases III and IV, during which the star slowly went back to its original 
appearance as the previously built disc dissipated). As discussed below, the timescales 
involved in the disc formation and dissipation are generally consistent with the timescales 
of viscous diffusion. 



2.7. Short-term variations 

Small-scale, short-term variations are quite common in Be stars, and possess a complex 
phenomenology (see Rivinius 2007, for a review). On one hand, some observed variations 
(e.g. line profile variations due to non-radial pulsations) are associated with the photo- 
sphere proper; others, on the other hand, are thought to originate from the very base of 
the disc and are, therefore, the manifestation of the physical process(es) that is (are) feed- 
ing the disc (e.g. short-term V/R variations of emission lines). To date, /x Cen remains 
the only system in which the ejection mechanism have been unambiguously identified 
(in this case non-radial pulsations, Rivinius et al. 1998). In the viscous diffusion theory 
outlined below, it is assumed that matter is ejected by the star and deposited at the 
inner boundary of the disc with Keplerian or super-Keplerian speeds; clearly, the current 
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Figure 2. Light curve (left panel) and CMD (right panel) of the variation of the star OGLE 
005209.92-731820.4 (de Wit 2006, reproduced with permission). 



status of the theory is still unsatisfactory inasmuch as the fundamental link between the 
disc and the photosphere proper (i.e., the feeding mechanism) is still unknown. 

Studies of short-term variations associated with the disc (e.g., Rivinius et al. 1998) are 
generally consistent with the following scenario: mass injection by some mechanism that 
causes a transient asymmetry in the inner disc that manifests itself by, e.g., short-term 
V/R variations, followed by circularization and dissipation in a viscous scale (see below). 
Here, again, viscosity seems to be playing a major role in shaping the evolution of ejecta. 



3. Viscous Decretion Disc Models 

The observational properties outlined above give strong hints as to the ingredients a 
successful theory for the structure of the Be discs must have. The theory must explain 
the keplerian rotation, slow outflow speeds, small geometrical thickness and must also 
account for the timescales of disc build-up and dissipation. 

The only theory to date that satisfies those requirements is the viscous decretion disc 
model (Lee et al. 1991). This model is essentially the same as that employed for pro- 
tostellar discs (Pringle 1981), the primary difference being that Be discs are outowing, 
while pre- main-sequence discs are inflowing. In this model, it is supposed that some 
yet unknown mechanism injects material at the Keplerian orbital speed into the base 
of the disc. Eddy /turbulent viscosity then transports angular momentum outward from 
the inner boundary of the disc (note that this requires a continual injection of angular 
momentum into the base of the disc). If the radial density gradient is steep enough, angu- 
lar momentum is added to the individual fluid elements and they slowly move outward. 
To critically test such decretion disc models of Be stars against observations, we must 
determine the structure of the disc from hydrodynamical considerations. 



3.1. Viscous discs fed by constant decretion rates 

There are many solutions available for the case of a viscous disc fed by a constant de- 
cretion rate, all agreing in their essentials. Below we outline the mains steps of one 
possible derivation (Bjorkman 1997) to the goal of discussing the properties of the 
solution and confronting them to the observations. More detailed analyses as well as 



different approaches to the problem can be found in Lee et al. (1991) |Okazaki (2001) 



Bjorkman k Carciofi (2005) Okazaki (2007) Jones et al. (2008) and Carciofi k Bjorkman (2008) 
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3.2. Hydrostatic Structure 

Our goal is to write and solve the Navier Stokes fluid equations in cylindrical coordinates 
(vj, <f>, z) which in the steady-state case have the following general form 

+ 1 + JL w = o , (3.i) 

dv^ dv w dv w w| l dP 

h —-^T+v z -; - = — TT- + , (3.2) 

ovj vj o<p oz vj p ovj 

dvj, w dvj, dv^ v^v,p 1 OP 

^-k^ + — -of + w * o H = "oT + /0 ! 3 - 3 

aro vj 0(f) oz vj pvj o<p 

dv z dv z dv z ldP 

+ — -5-r = — "5-+A 5 ( 3 - 4 ) 

aro vj 0(f> oz p oz 

where p is the gas mass density, P is the pressure and / w , and / 2 are the components 
of the external forces acting on the gas. 

Let us initially ignore any viscous effects (inviscid disc) and assume that the only force 
acting on the gas is gravity. If we further assume circular orbits — u ro = 0, 7^ 0, and 
v z = 0, an assumption that will be droped later on when viscosity is included — the only 
non-trivial fluid equations are the vj- and z-momentum equations (Eqs. 13.21 and I3.4j) . 
which take the form 

(3.5) 

p ovj vj 

~~tr = fz , (3.6) 

p oz 

where the external force components are given by the gravity of the spherical central star 

GMw 

Jro (vj 2 + z 2 ) 3 / 2 ' [ ' 

f> = - {m ? +z 2 )V 2 ■ ( 3 - 8 ) 

To specify the pressure, we introduce the equation of state, P = c 2 p, where c s = 
(/cT) 1 / 2 (/imH) -1 ^ 2 . In this last expression, k is the Boltzmann constant, p is the gas 
molecular weight and toh is the mass of the hydrogen atom. 

In the thin-disc limit (z -C vo), we obtain 



= Krit (RM 1/2 , (3.9) 

(3.10) 



dln(c 2 p) _ V 2 it Rz 



dz c 2 vj 3 

The above equations mean that the disc rotates at the Keplerian orbital speed and is 
hydrostatically supported in the vertical direction. To determine the vertical disc struc- 
ture we must solve Eq. (|3.10l) . Assuming an isothermal disc one obtains 

p(vj, z) = p (vj) exp [-0.5(z/Hf] , (3.11) 

where po the disc density at the mid-plane (z = 0), and the disc scale height is given by 

H(vj) = (cs/v^vj. (3.12) 

Since v<f, oc tn~ 5 , we obtain the familiar result that for an isothermal disc the scaleheight 
grows with distance from the star as H oc vj 1,5 . 
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As we shall see below, it is useful to express p(w, z) in terms of the disc surface density, 
£, written as 



p{vj, z) dz — V^nHpo . (3.13) 

-oo 

p(w, z) = I**) exp [-0.5(z/Hf] . (3.14) 



Thus 



3.3. Viscous Outflow 

Clearly, in Eq. (|3.14j) the disc density scale, J^(ur), is completely undetermined, because 
for a inviscid Keplerian disc we can choose to put an arbitrary amount of material at a 
given radius. To set the density scale, we must include a mechanism — viscous diffusion 
— to transport material from the star outwards. 
Viscous flows grow in a viscous diffusion timescale 

T diS = m 2 /v, (3.15) 

where v is the kinematic viscosity. One problem that has already been noted long ago 
(Shakura & Sunyaev 1973) is that for molecular viscosity the diffusion timescale is much 



too long. Shakura & Sunyaev (1973) appealed instead to the so-called eddy (or turbulent) 
viscosity, which they parameterized as 

v = ac s H , (3.16) 

where < a < 1 . The a parameter describes the ratio of the product of the turbulent 
eddy size and speed to the product of disc scaleheight and sound speed. In other words, it 
is assumed that the largest eddies can be at most about the size of the disc scaleheight and 
that the turnover velocity of the eddies cannot be arger than the sound speed (otherwise, 
the turbulence would be supersonic and the eddies would fragment into a series of shocks) . 
With this value of the viscosity, the viscous diffusion timescale becomes 

Tdiff = ^^VwR (3-17) 




(3.18) 

Studies of the formation and dissipation of the discs around Be stars find typical time 
scales of months to a few years (e.g., Wisniewski et al. 2010); therefore, an a of the order 
of 0.1 or larger is required to match the observed timescales (see Sect. HJ. 

If we add viscosity to the fluid equations, we can still assume that the disc is axisym- 
metric and that the vertical structure is hydrostatic (v z = 0). However, the presence of 
an outflow implies that w ro ^ 0. The w- and z- momentum equations are the same as 
before, so v$ and p are the same as in the pure Keplerian case [eqs. (|3.9j) and (|3 . 14[) ] . 

Two fluid equations remain to be solved, the continuity equation, Eq. (|3.1[) . and the </>- 
momentum equation, Eq. (|3.3p . The continuity equation, written in terms of the surface 
density, 

d 

— (27ra7&; ro ) = (3.19) 
ow 

means that the mass decretion rate, M = 2tivjYiV^, is a constant (independent of w). 
The viscous outflow speed is given by 
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This (^momentum equation now is more complicated because viscosity exerts a torque, 
which is described by the viscous shear stress tensor, 7r ro 0. Including this shear stress, 
the 0-momentum equation becomes 

dv^ v^ V(j> 1 2 

^U7-^- £ - H = T^— (^ 7T ro 0) , (3.21) 

aro zu pm* otxj 

where 

d{v^/w) 3 2 

7T ro = — = -^ aC sP ■ ( 3 - 22 ) 

Multiplying Eq. (|3.21|) by pw 2 and integrating over <j> and z, we find 

M-^i™*) = ^(T) , (3.23) 

where 

/>oo 

T = I WK^^-Kwdz = —2>-Kac 2 s w 2 Yi (3.24) 



is the viscous torque. The <p- momentum equation, Eq. p.21| . expresses the fact that the 
change in the angular momentum flux — vjv^ being the specific angular momentum — 
is given by the gradient of the viscous torque. Since the continuity equation implies that 
M is constant, we integrate Eq. (|3.21l) over w to obtain 

Tim) = MV clit VvjR + constant . (3.25) 
Substituting Eq. (|3.24p into Eq. (|3.25p and solving for the surface density we find 

SH = ^ (^) [(^M 1/2 l] • (3-26) 

Eq. (|3.26p describes the surface density of an unbounded disc (i.e. a disc that is allowed 
to grow indefinitely). The integration constant Ro is a parameter that depends on the 
integration constant of Eq. (|3.25|) and is related to the physical size of the disc; for time- 
dependent models, such as those of Okazaki (2007) [ Rq grows with time and thus Ro is 
related with the age of the disc. 

3.4. Properties of the Solution 

We have now completed the hydrodynamic description of an isothermal, unbounded 
viscous decretion disc. To fully determine the problem one must specify the decretion rate, 
M, the value of a, and the disc age (or size). Assuming that the disc is sufficiently old, 
Ro ^> R, in which case Eq. (|3.26[) becomes a simple power-law with radius, E(tJ7) cx vo~ 2 . 
From Eqs. (|3.14[) and (|3.12[) we obtain that the isothermal disc density profile is quite 
steep, p oc vu~ 3 ' 5 . 

Another important property of isothermal viscous discs can be readily derived from 
Eq. (|3.12[) . Since is much larger than the sound speed (the former is of the order of 
several hundreds of km/s whereas the later is a few tens of km/s), the disc scaleheight 
is small compared to the stellar radius, i.e., the disc is geometrically thin. Finally, from 
Eq. p. 201) we find that, for large discs, the radial velocity is a linear function of the radial 
distance, v ro cx w. 

3.5. Temperature structure 



The temperature structure of viscous discs were investigated by Carciofi & Bjorkman (2006) 



and Sigut & Jones (2007) They found that those discs are highly nonisothermal, mainly 
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Figure 3. Top: Vertically averaged disc temperature. Bottom: Local power-law index of the 
density profile. Left: Results for a low-density model (po = 4.2 x 1 0~ 12 g cm" 3 ). Right: Results 
for a high-density model (po = 3.0 x 10 -11 g cm -3 ). Adapted from |Carciofi fc Bjorkman~ (2008) 



on their denser inner parts. An example of the temperature structure is shown in Fig. [31 
The temperature initially drops very quickly close to the stellar photosphere; when the 
disc becomes optically thin vertically the temperature rises back to the optically thin ra- 
diative equilibrium temperature, which is approximately constant, as in the winds of hot 
stars. The density is the most important factor that controls the temperature structure: 
for high density models the amplitude of the temperature variations is much larger and 
the nonisothermal region extends much farther out into the disc. 



3.6. Non-isothermal effects on the disc structure 

From Eqs. (13.141) and (|3.26p we see that both the viscous diffusion and the vertical hy- 
drostatic solutions depend on the gas temperature; one can expect, therefore, that the 
complex temperature structure of the disc might have effects on the disc density struc- 
ture. Carciofi & Bjorkman (2008) and Sigut et al. (2009) calculated the vertical density 



structure of a disc in consistent vertical hydrostatic equilibrium. They found that the 
temperature decrease causes the disc to collapse, becoming much thinner in the inner 
regions. This collapse redistributes the disc material toward the equator, increasing the 
midplane density by a factor of up to 3 relative to an equivalent isothermal model. 



Carciofi & Bjorkman (2008) investigated, in addition, how the temperature affects the 



viscous diffusion. The combination of the radial temperature structure, disc scaleheight, 
and viscous transport produces a complex radial dependence for the disc density that 
departs very much from the simple n = 3.5 power-law. As shown in Fig.[3j the equivalent 
radial density exponent varies between n = 2 in the inner disc to n = 5 near the temper- 
ature minimum, eventually rising back to the isothermal value n = 3.5 in the outer disc. 
We conclude that non-isothermal effects on the viscous diffusion may account for the at 
least part of the large scatter of the index n reported in the literature (Sect. 12. 3|) . 
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Figure 4. Left: Fit to the SED of £ Tau. The dark grey lines are the observations and the black 
lines the model results. The light grey line corresponds to the unatten uated stellar SED and 
gives a measure of how the disc affects the emergent flux. Adapted from |Carciofi et al. (2009) | 



3.7. Two test cases: C Tau and \ Oph 

The model described above makes several predictions about the disc structure that are 
in qualitative agreement with the observations (Sect. [5]): viscous discs are geometrically 
thin, rotate in a near keplerian fashion (v^ oc vj^ 1 / 2 and v ro <C v^) and have a steep 
density fall-off (n = —3.5 in the isothermal case). This model and its predictions must 
now be quantitatively compared to observations. 

A successful verification of the viscous decretion disc model have been obtained in the 
case of the Be star ( Tau (Carciofi et al. 2009). This star is particularly suitable for such 
study because it had shown little or no secular evolution in the past 18 years or so (Stefl 
et al. 2009); therefore, a constant mass decretion rate, as assumed above when deriving 
the disc structure, is a good approximation for this system. Some of the results obtained 
for £ Tau are shown in Fig. [4] The radial dependence of the disc density, temperature, 
and opening angle all affect the slope of the visible and IR SED (Waters 1986, Carciofi 
& Bjorkman 2006), as well as the shape of the intrinsic polarization. Therefore, the fact 
that the model reproduces the detailed shapes of SED and spectropolarimetry represents 
a non-trivial test of the Keplerian decretion disc model. 

Using a somewhat different model, Tycner et al. (2008) successfully fitted several ob- 
servations of the Be star \ Oph, including high angular resolution interferometry. In their 
modeling the index n of the density power-law is a free parameter. Interestingly, their 
best fitting model was for a much flatter density law (n = 2.5). Whether this flatter den- 
sity profile can be accommodated by including some other physical process in addition 
to viscosity (e.g., tidal effects by a binary) or is an indication that the viscous decretion 
disc model is not a good model for this system remains to be verified by further analysis. 



4. Viscous Discs Fed by Non-constant Decretion Rates 

All the discussion so far have been focused on the problem of a disc fed by a constant 
decretion rate that grows to a given size in a viscous timescale. Clearly, those models can 
only be applied to objects, such as £ Tau, that went through a sufficiently long and stable 
decretion phase. A different approach is needed if one wants to investigate dynamically 
active systems such as the one shown in Fig. [2] 

Okazaki (2007) | Haubois et al. (2010) and Jones et al. (2008) described the solution of 



the viscous diffusion problem for systems with non-constant mass decretion rates. Fig. [5] 
shows a series of models of the lightcurve of the Be star 28 CMa (Carciofi et al., in 
prep.), that uses the computer code SINGLEBE by Atsuo Okazaki to solve the ID viscous 
diffusion problem and the hdust code to calculate the emergent spectrum (Haubois et 
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Figure 5. The visual light curve of 28 CMa fitted with a dynamical viscous decretion disc 
model with different values of the viscosity parameter, a (Stefl et al. 2010) 

al. 2010). The simulation begins after a few years of disc evolution to account for the 
previous disc build-up. At t = 2003.6 mass decretion is turned off and the system evolves 
passively until t = 2008.8, when the recent outburst started (Stefl ct al. 2010). The model 
with a — 0.9 is the one that reproduces best the lightcurve at all phases. This is, to our 
knowledge, the first time the viscosity parameter is determined for a Be star disc. 

5. conclusions 

This paper discussed the basic hydrodynamics that determine the structure of viscous 
decretion disc models. Those are, to date, the most satisfactory models for Be discs 
because they can account quantitatively for most of their observational properties, namely 
the fact that they are geometrically thin, have Keplerian rotation and steep radial density 
profiles. It has been recently shown that the viscous decretion disc model can also explain 
the temporal evolution of dynamically active systems such as 28 CMa. Studies of such 
system allows for the determination of the viscosity parameter a and the mass decretion 
rate, M, two quantities that are very difficult to determine otherwise. 
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Discussion 

PULS: Excellent talk! How are the occupation numbers required for the profile synthesis 
calculated? 

Carciofi: The computer code hdust calculates the occupation numbers by solving the 
coupled problem of the radiative transfer, radiative equilibrium and statistical equilib- 
rium in the NLTE regime. Because the code uses the Monte Carlo method, it can handle 
arbitrary 3D density and velocity distributions. 

Lee: What is the timescale of the V/R variability? 

Carciofi: If the question refers to cyclic V/R variability, such as observed in £ Tau, the 
timescale is of a few years 



